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Abstract. A set of 45 dust envelopes of carbon stars 
has been modeled. Among them, 34 were selected accord- 
ing to their dust envelope class (as suggested by Sloan, 
Little-Marenin & Price, 1998) and 11 are extreme car- 
bon stars. The models were performed using a code that 
describes the radiative transfer in dust envelopes con- 
sidering core/mantle grains composed by an a-SiC core 
and an amorphous carbon (A.C.) mantle. In addition, we 
have also computed models with a code that considers 
two kinds of grains - a-SiC and A.C. - simultaneously. 
Core-mantle grains seem to fit dust envelopes of evolved 
carbon stars, while two homogeneous grains are more able 
to reproduce thinner dust envelopes. Our results suggest 
that there exists an evolution of dust grains in the carbon 
star sequence. In the beginning of the sequence, grains are 
mainly composed of SiC and amorphous carbon; with dust 
envelope evolution, carbon grains are coated in SiC. This 
phenomena could perhaps explain the small quantity of 
SiC grains observed in the interstellar medium. However, 
in this work we consider only a-SiC grains, and the in- 
clusion of /3-SiC grains can perhaps change some of there 
results. 

Key words: Stars:carbon - circumstellar matter - radia- 
tive transfer - SiC grain 



1. Introduction 

Asymptotic giant branch (AGB) stars are often sur- 
rounded by circumstellar dust shells. The chemical com- 
position of these media reflects that of the stellar photo- 
sphere. Thus, carbon-rich grains, such as amorphous car- 
bon (A.C), are one of the expected components of the 
dust envelopes around carbon stars. In addition, almost all 
these stars show an emission feature around 11.3/im due 
to Silicon carbide (SiC) grains which are also condensed 
there. The existence of SiC grains in the atmospheres of 
carbon stars was predicted at first on the basis of chem- 
ical equilibrium calculations by Gilman (1969). This pre- 
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diction was supported by the observations of Hackwell 
(1972), Treffers & Cohen (1974) and Forrest et al. (1975). 
Nowadays infrared satellites provide several characteristic 
features of these grains. 

After the observations with IRAS satellite, several 
works have dealt with the classification of dust envelopes 
around carbon stars. Little-Marenin et al. (1987) have dis- 
covered 176 new carbon stars using the feature at 11.3/im 
as a selection factor. Willems (1987) has analyzed 304 such 
objects, suggesting a SiC-index to the stars which present 
the emission feature varying between 11.2-11.6/im. Pa- 
poular (1988) has classified carbon dust envelopes using 
the 11.3/im feature and a secondary feature at 8.6/im. The 
sample of Chan & Kwok (1990) was composed of 356 ob- 
jects which were classified in two distinct classes. Accord- 
ing to the authors, the difference between the classes is 
due to an evolution of a-SiC and /3-SiC particles. More 
recently Sloan, Little-Marenin & Price (1998, hereafter 
SLMP) proposed a classification of 89 carbon-rich stars 
in 6 different types based on their infrared emissions. All 
stars of their sample show the SiC feature at 11.3/im. 

In general, differences between each class can be inter- 
preted as an evolution of the dust envelope itself due to an 
increasing amount of grains; consequently, optical depth 
also increases, affecting radiative transfer and the emis- 
sion feature. On the other hand, the resonance features 
of small SiC grains are very sensitive to size, morphology 
and chemical composition of impurities in the surrounding 
medium (Bohen & Huffman 1998). This fact suggests that 
the variety of emission features assigned to SiC grains is 
probably related to the formation process of the dust in 
circumstellar envelopes, reflecting the physical and chem- 
ical conditions. The classification criteria cited above are 
useful, but if one is interested in a deeper insight into the 
nature of the circumstellar dust grains, it is necessary to 
solve the radiative transfer problem, and to reproduce the 
features seen in the mid-IR LRS together with the overall 
behavior of the spectral energy distribution. 

Several authors have calculated the radiative transfer 
in circumstellar dust shells (CDS) of carbon stars (Chan & 
Kwok 1990; Lorenz-Martins & Lefevre 1993, 1994; Groe- 
newegen 1995; Bagnulo, 1996; Bagnulo, Doyle and An- 
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dretta 1998). Lorenz-Martins & Lefevre (1994) have em- 
ployed the Monte Carlo method for solving the radiative 
transfer for two species simultaneously. SiC grains were 
supposed to form closer to the star than graphite grains 
(McCabe 1982). Correlations between SiC/A.C. ratio and 
extinction opacity as well as SiC/A.C. ratio and period 
of luminosity were found. These correlations indicate that 
the quantity of SiC grains relative to amorphous carbon 
grains decreases with carbon star evolution. Mass loss 
from cool stars is the major source of refractory grains in 
the interstellar medium. Carbon stars provide carboneous 
material and extreme carbon stars, with their large mass- 
loss rate, are the main contributors. A study by Whittet 
et al. (1990) shows that there is a low fraction (< 5%) 
of Si in the form of SiC in the interstellar medium. The 
enrichment rate for a kind of grain depends on the com- 
position of the dust envelopes of stars which have high 
mass-loss rate. Lorenz-Martins & Lefevre (1994) also sug- 
gested that SiC grains are the minor component in carbon 
star envelopes. In addition, stars which have low SiC/A.C. 
ratios (0.01-0.06) have also high mass-loss rates (1.5.10~ 4 
- IRC+10216 — 6.1.10" 6 M Q /yr - IRAS 15194-5115). 
These results can perhaps explain the small quantity of 
SiC grains in the interstellar medium. Alternatively, the 
weakness of the SiC absorption could be due to the fact 
that SiC particles are embedded in thick carbon mantles. 
In fact, Kozasa et al. (1996) have shown that the nucle- 
ation of SiC precedes that of carbon grains and may lead 
to the formation of dust grains consisting of a SiC core 
and a carbon mantle. Moreover they proposed that such 
core/mantle grains are the most reasonable candidates to 
reproduce the feature seen around 11.3/im. They also sug- 
gested that radiative transfer calculations should be used 
in order to verify this model. 

The main purpose of the present work is to present 
dust envelope models considering core/mantle grains com- 
posed by an a-SiC core and an A.C. mantle. We also 
compare the data with models consisting of a-SiC and 
A.C. homogeneous grains. Two samples of stars have been 
considered. The first one contains objects classified by 
SLMP. They were analyzed aiming to verify the classi- 
fication proposed by the authors. The second sample con- 
tains 11 extreme carbon stars, which have thicker dust 
envelopes and a higher mass loss rate. The utilization of 
core/mantle grains is satisfactory to describe some stars 
but most of them are better reproduced using a two ho- 
mogeneous grains model. 

2. Core/mantle grains and the method 

Besides carbon grains of different structures (e.g., 
graphite, amorphous carbon), SiC is the most important 
species for late-type carbon stars. SiC is one of the most re- 
fractory materials that may condense under the conditions 
of a carbon- rich chemistry. In such stars, almost all oxy- 
gen is chemically blocked in the CO molecule and, among 



the more abundant chemically active elements, carbon and 
SiC are the possible condensates. An important problem 
in the formation of SiC under the conditions present in 
circumstellar shells is that the abundance of SiC molecule 
is quite low. Beck (1992) has considered non-equilibrium 
effects and has shown that solid SiC can be stable against 
evaporation at temperatures below 1400K. He also sug- 
gested that SiC may form on the surface of preexisting 
carbon grains instead of being the primary condensate 
at very high temperatures. On the other hand, McCabe 
(1982) has shown that SiC particles can be formed at high 
temperatures due to the greenhouse effect. Following this 
suggestion, Kozasa et al. (1996) have demonstrated that 
the nucleation of SiC grains always precedes that of car- 
bon grains when the non-LTE effect, i.e., the difference be- 
tween the temperature of gas and small clusters, is taken 
into account. 

We have considered this latest suggestion and modeled 
45 envelopes of carbon stars using core/mantle grains con- 
sisting of a a-SiC core and a A.C. mantle. The method cm- 
ployed is an improved version of that described by Lorenz- 
Martins & Araiijo (1997). We have modified the code to 
include a new option about grain properties. The absorp- 
tion and scattering efficiencies, as well as the albedo, were 
calculated using the Mie theory for core/mantle grains 
(e.g. Bohren & Huffman, 1984; Hoyle & Wickramasinghc, 
1991) and optical constants (or dielectric functions) tab- 
ulated in the literature. The optical constants which we 
have used are the ones determined by Pegourie (1988) for 
a-SiC, and by Rouleau & Martin (1991) for amorphous 
carbon. 

The propagation of stellar and grain radiative energy 
is simulated photon by photon following a Monte Carlo 
scheme. For each interaction between a "photon" and a 
grain, a fraction of the energy is stored (absorption) and 
the remaining part is scattered according to the scatter- 
ing diagram. The stellar radiation leads to an initial dis- 
tribution of dust temperature and the thermal radiation 
from grains is simulated, giving after several iterations the 
equilibrium temperature. Computations give the spectral 
repartition of the total flux and of its different components 
(direct, scattered, emitted), and the temperature law for 
the grains. For more details see Lorenz-Martins & Lefevre 
(1993, 1994) and Lorenz-Martins & Araujo (1997). 

3. The sample 

SLMP have studied 89 carbon-rich stars and organized 
the dust emission in several classes. Red class contains 
only 3 stars which present a 11.3/im feature and a strong 
dust continuum. SiC class is the most numerous, with 40 
objects. They have the 11.3/im feature and a weak dust 
continuum. In the SiC+ they put 32 stars which show 
two features: the 11.3^m one and a weak feature at 8- 
9/im. These objects also have a weak dust continuum. The 
SiC++ class contains 6 stars with comparable 11.3/xm and 
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Table 2. Sample with evolved stars 



TD AQ 
1-TurVD 


^1 q in a 




°p- iyp 6 






05377+1346 


AFGL 


799 


C8,4 


372 c 


1,4 


05405+3240 


AFGL 


809 


C 


780 6 


5 


06012+0726 


AFGL 


865 


? 


696 c 


6 


06291+4319 


AFGL 


954 


c 




5 


06342+0328 


AFGL 


971 


c 


653 c 


6 


07098-2112 


AFGL 


1085 


N 


725 c 


6 


08088-3243 


AFGL 


1235 


c 


571 c 


2 


15082-4808 


AFGL 


4211 


? 




3 


19594+4047 


AFGL 


2494 


c 


783 6 


5 


20570+2714 


AFGL 


2686 


C8,5 


750 6 


5 


23257+1038 


AFGL 


3099 


C 


484 c 


5 



(1) Cohen (1984); (2) Epchtein et al. (1990); (3) Fouque et al. 
(1992); (4) Gehrz & Hackwell (1976); (5) Jones et al. (1990); 
(6) Le Bertre (1992); Periods: (o) Kholopov et al. (1987), (b) 
Jones et al. (1990), (c) Le Bertre (1992). 



8-9/im features. The five stars in the Broad 1 class have 
an unusual 11.3/xm feature profile with short- wavelength 
excess. Finally, the Broad 2 class contains only 3 stars 
which present an unusual 11.3/xm feature profile with long- 
wavelength excess. 

In Table 1 we present our sample of 34 stars taken 
from SLMP. We restrict our study to the objects with 
IR fluxes published in the literature. Fortunately, we have 
been able to obtain data from stars belonging to all dif- 
ferent classes. Table 1 lists the IRAS number (column 1) 
followed by the usual name (column 2). Spectral type and 
variability class are listed in columns 3 and 4 respectively. 
Column 5 gives the period and column 6 shows the enve- 
lope class attributed by SLMP. Some stars have a doubtful 
classification; they are designated by a colon. Finally the 
photometry used in the fit of the models is presented in 
the last column. Table 2 list the extreme carbon stars sam- 
ple. The columns are analogous to those of Table 1. The 
stars analyzed are variables and self-consistent models for 
them require data from similar phase of luminosity. How- 
ever this is difficult due to the scarcity of observations. We 
have worked with the photometry available in the litera- 
ture and considered the phase whenever possible. In order 
to minimize the uncertainties, we have used average LRS 
(IRAS, 1986) spectra for SiC emission. The 12, 25, 60, 
100 fim fluxes were taken from Gezari et al. (1987) and 
SIMBAD. 

4. Results 

Table 3 and Table 4 present the results of the best mod- 
els to our first and second samples, respectively. These re- 
sults were obtained considering a core/mantle grain. IRAS 
number (column 1) is followed by the temperature of the 
central star (Teff in K) in the second column. Third and 
fourth columns present inner (R x in R*) and outer (R 2 



in R*) envelope radii, respectively. The dimension of the 
a-SiC core (csic in A) of the grains is given in the fifth col- 
umn, followed by the dimension of the amorphous carbon 
mantle (ula.c. m A) in the sixth column. Finally, the opti- 
cal depth (t) at l[im and the abundance ratio between a- 
SiC and amorphous carbon grains are given, respectively, 
in columns 7 and 8. In Table 3 we have added one last 
column with the SLMP's envelope classes. 

In addition, we have calculated models for the stars 
in the sample using a code with two homogeneous grains. 
Some stars of this sample were analyzed in Lorenz-Martins 
& Lefevre (1994), where the authors describe the method. 
The differences between parameters of both codes are the 
size of the grains and the way in which the SiC/A.C. ratios 
were calculated. In the core/mantle models, we consider 
the core (csic) and mantle (mA.c.) size. SiC/A.C. ratios 
in this method are obtained by mass, based on the man- 
tle and core size, and we find the corresponding value as 
obtained in the two homogeneous grains model. Others 
parameters are obtained the same way as in the two ho- 
mogeneous grains code. 

The features seen in the mid-IR LRS must be repro- 
duced and the properties of complete CDS must be deter- 
mined simultaneously. In order to fit the dust emission, 
we have calculated grids of about fifty models for each 
star, and we inspect visually the model which best repro- 
duces the complete CDS. We pay special attention to the 
LRS feature and accept errors of about 15 per cent in all 
parameters (such as grain size, effective temperature, op- 
tical depth...) to fit this emission feature; in fact it is this 
feature that defines the best model. 

4.1. Results for SLMP stars 
SiC Class 

We analyzed 9 out of 40 SLMP objects in this class. 
We found that effective temperatures vary between 2200 
K and 2700 K. Inner radii vary between 3.7R* and 5. OR*, 
while most of the outer radii are 1000R*. It must be kept 
in mind that the results are not very sensitive to this last 
parameter, as has been pointed out in previous works. The 
a-SiC core (csic) grains vary between 280 to 550 A, and 
the amorphous carbon mantle (mA.c.) between 700 and 
1050 A. Optical depths for this class vary between 0.10 
and 2.5, with two exceptions: U Hya (r = 0.03) and R 
For (t = 7.0). Finally SiC/A.C. ratios vary from 0.06 to 
0.20. In Figures l(a)-(d) we show out fits to RV Aqr and 
Y Tau. Solid lines represent the core/mantle grain model 
and dashed lines the two homogeneous grains model. Fig- 
ures 1 (b) and 1 (d) show an enlarged view centered on the 
11.3/zm feature. In all cases, best fits to this class were 
obtained using the two homogeneous grains model. 

SiC+ Class 

We have modeled 7 out 32 SLMP SiC+ stars. Accord- 
ing to our results, the temperatures of central stars vary 
from 2400 K to 2650 K. Outer dust envelope radii are the 



4 S.Lorenz-Martins et al.: Modeling of C stars with core/mantle grains 



Table 1. Sample of SLMP's stars 





Name 


op-iype 


var. 


Period 


ii/iiv. Class 




00172+4425 


VX And 


C4,5J 


SRa 


369 a 


SiC++ 


9 


01246-3248 


R Scl 


C6,5 


SRb 


370 a 


SiC++ 


1,7 


02270-2619 


R For 


C4,3e 


Mira 


388 c 


SiC 


8 


03075+5742 


C* 131 


C4,5J 


Lb 


— 


N: 


9 


03374+6229 


U Cam 


C6,4 


SRb 


— 


SiC+ 


9 


04459+6804 


ST Cam 


C5,4 


SRb 


300 a 


SiC+: 


9 


04483+2826 


TT Tau 


C7,4 


SRb 


166 a 


SiC+: 


9 


04573-1452 


R Lep 


C7,4e 


Mira 


427 a 


SiC 


2,8 


05028+0106 


W Ori 


C5,4 


SRb 


212 a 


SiC+ 


2,9 


05418-4628 


W Pic 


C 


Lb 


— 


SiC++: 


7 


05426+2040 


Y Tau 


C6,4 


SRb 


241 a 


SiC 


7 


05576+3940 


AZ Aur 


C8 


Mira 


416 a 


Br2 


9 


06225+1445 


BL Ori 


C6,3 


Lb 


— 


Br2 


3,9 


06331+3829 


UU Aur 


C7,4 


SRb 


23a 


SiC 


2,9 


06529+0626 


CL Mon 


C6,3e 


Mira 


497" 


SiC 


9 


07045-0728 


RY Mon 


C5,5 


SRa 


278" 


SiC+ 


9,12 


07057-1150 


W CMa 


C6,3 


Lb 


— 


SiC+: 


7,9 


07065-7256 


R Vol 


Ce 


Mira 


454" 


SiC 


5 


08538+2002 


T Cnc 


C5,5 


SRb 


482" 


SiC++ 


2,3 


09452+1330 


IRC+10216 


C9,4 


Mira 


649 c 


Red 


11 


10329-3918 


U Ant 


C5,3 


Lb 


— 


SiC+: 


7 


10350-1307 


U Hya 


C5,3 


Lb 




SiC 


2,5,9 


10491-2059 


V Hya 


C6,5 


SRa 


531 a 


Red 


2,5,9 


12226+0102 


SS Vir 


C6,3 


SRa 


364" 


Brl 


9 


12427+4542 


Y CVn 


C5,5J 


SRb 


158" 


SiC+: 


2,9,11 


12447+0425 


RU Vir 


C8,le 


Mira 


433" 


SiC 


5 


12544+6615 


RY Dra 


C4,5J 


SRb 


200 a 


SiC+ 


2,9 


15094-6953 


X TrA 


C5,5 


Lb 




SiC+ 


5,11,12 


18306+3657 


T Lyr 


C6,5J 


Lb 




SiC++ 


9 


19017-0545 


V Aql 


C6,4 


SRb 


353 a 


SiC+ 


5,9 


19555+4407 


AX Cyg 


C4,5 


Lb 




SiC+: 


9 


21032-0024 


RV Aqr 


C6,3e 


Mira 


454" 


SiC 


5 


21399+3516 


V460Cyg 


C6,4 


SRb 


263 a 


SiC+: 


2,9 


23587+6004 


WZ Cas 


C9,2J 


SRb 


186" 


Brl 


3,10 



(1) Bagnulo et al. (1998) (2) Bergeat et al. (1976); (3) Bergeat & Lunel (1980); (4) Epchtein et al. (1990); (5) Fouque et al. 
(1992); (6) Jones et al. (1990); (7) Kerschbaum, Lazaro& Habison (1996); (8) Le Bertre (1992); (9) Noguchi et al.(1981); (10) 
Noguchi & Akiba (1986); (11) Nyman et al. (1992); (12) Walker (1976). Periodos: (a) Kholopov et al. (1987), (b) Jones et al. 
(1990), (c) Le Bertre (1992). (i) SLMP 



same for all stars (R2 = 1000 R*) and inner radii vary 
between 3.7 and 6.5 R*. Sizes of mantle grains (m^.c.) 
show a great dispersion: 400 A to 1200 A. The same occurs 
with the sizes of a-SiC core (cgic) which vary from 155 A 
to 500 A. Optical depths have values between 0.03 and 
0.50, and SiC/A.C. ratios between 0.06 and 0.26. Almost 
all stars were best described using the two homogeneous 
grains model; the unique exception is V Aql, which is bet- 
ter reproduced with a core/mantle grain code. Figure 2 
shows best fits to RY Mon and U Cam. 

As a rule, we can say that the envelope of SiC and 
SiC+ stars are nicely described by the existence of a-SiC 
and A.C. homogeneous grains simultaneously. 



We analyzed 5 out of 6 SLMP objects. According to 
our models, effective temperatures are either 2200K or 
2400K. Almost all inner and outer radii have the same 
values: Ri = 5R* and R2 = 1000R*. Mantle sizes (niA.c.) 
vary between 700 and 1200 A, and core sizes (csic) be- 
tween 300 and 500 A. Optical depths vary between 0.02 
and 0.20, which indicates very thin dust envelopes, and 
SiC/A.C. ratios between 0.04 and 0.10. Due to the absence 
of optical constants to describe the 8-9 /im and 13 /zm 
features, it is difficult to choose between the core/mantle 
grain or the two homogeneous grains model. We will dis- 
cuss these results in the next section. Figure 3 shows two 
stars of this class. 



SiC++ Class 



Broad 1 and Broad 2 Classes 
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Wavelength (**m) Wavelength (urn) Wavelength Um) Wavelength (nm) 



Fig. 1. Best models for SiC class. In Figure l(a)-(d) we 
plotted the best models using core/mantle grains (solid 
line) and two homogeneous grain model (dashed line). Tri- 
angles represent the photometric data. In Figure 1(b) and 
1(d) we show an enlarged view of the 11.3 /zm emission. 



10" s | 1 i i i m i 1 r i i n il, i—i | 1 i i i u i h 10" 1 




Wavelength (**m) Wavelength (urn) 



Fig. 2. This figure shows the best models for RY Mon 
and U Cam, which are in the SiC+ class. In Figure 2(a)- 
(d) we plotted models using core/mantle grain (solid line) 
and two homogeneous grain model (dashed line) . Triangles 
represent the photometric data. In Figure 2(b) and 2(d) 
we show an enlarged view of the 11.3 ^m emission. 



Fig. 3. This figure shows best models for TCnc and R Scl, 
which arc in SiC++ class. In Figure 3(a)-(d) we plotted 
models using core/mantle grain (solid line) and two ho- 
mogeneous grain model (dashed line). Triangles represent 
the photometric data. In Figure 3(b) and 3(d) we show an 
enlarged view of the 11.3 fxm emission. 

We have modeled 2 out 5 SLMP Broad 1, and 2 out 
3 SLMP Broad 2 stars. Our results indicate very simi- 
lar properties for both classes. The temperatures of the 
central stars are in the range of 2200K to 2700K. Core 
and mantle sizes vary between 400-500 A and 850-1200 A, 
respectively. They have thin envelopes with analogous di- 
mensions. SiC/A.C. ratios varying between 0.07 and 0.20 
were found. We have obtained our best fits using a two 
homogeneous grains code for 3 stars. WZ Cas presents an 
absorption at about 14 //m, which is usually seen in J-type 
carbon stars, and it's difficult to distinguish between the 
models. SS Vir shows features at about 8 //m and 14 /zm, 
too. Figure 4 presents best fit models to SS Vir {Broad 1 
class) and AZ Aur {Broad 2 class). 

Red Class 

We have modeled 2 out of 3 SLMP stars: IRC+10216 
and V Hya. Both objects are very well studied and sup- 
posed to have asymmetrical dust envelopes. They are be- 
lieved to be in the latest stages of stellar carbon evolution. 
We will discuss them below. 

(a) IRC+10216 

IRC+10216 has been extensively observed at optical, 
radio and infrared wavelengths. The central object is a 
long-period variable, with a period of 640 days, and it is 
commonly considered to be a late-type carbon star. Its 
envelope has been continuously surveyed and 380 molec- 
ular lines were detected, of which 317 have been identi- 
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Fig. 4. This figure shows best models for SS Vir, which 
belongs to Broad 1 class and AZ Aur to Broad 2 class. In 
Figure 4(a)-(d) we plotted models using core/mantle grain 
(solid line) and two homogeneous grain model (dashed 
line). Triangles represent the photometric data. In Figure 
4(b) and 4(d) we show an enlarged view of the 11.3 ^m 
emission. 



fied (Cernicharo, Guelin& Kahane 2000). Deep B and V 
image-bands, reveal its extended circumstellar envelope 
in the dust scattered and show an episodic mass loss rate. 
The circumstellar envelope is roughly spherically symmet- 
rical but it is likely to be composed of a series of discrete, 
incomplete, concentric shells (Mauron & Huggins 1999). 

IRC+10216 was modeled by several authors. Michell 
& Robinson (1980) have used graphite grains and radia- 
tive transfer calculations were performed considering, as 
usual, a spherically symmetric envelope with a central 
star. Rowan-Robinson & Harris (1983), Lc Bcrtrc (1987, 
1988b), Martin & Rogers (1987) and Griffin (1990) treated 
radiative transfer in the envelope of this star in a simi- 
lar way. Lorenz-Martins & Lefevre (1993,1994) have mod- 
eled this star by considering a two homogeneous grain 
model consisting of a-SiC and amorphous carbon grains 
simultaneously, as already cited. More recently, Groenewe- 
gen (1997) has computed a spherically symmetrical dust 
model and suggested that IRC+10216 is in the latest 
phases of carbon stars evolution, like V Hya. 

The effective temperature of IRC+10216 is 2100 K 
and our results lead to an extensive (Ri=5.6R» and R2 
= 8000R*) but thick (r = 10) dust envelope. The size of 
the amorphous carbon mantle grains (nt rm A.c.) was 500 A 
with an a-SiC core of 110 A. Our best model was obtained 



considering core/mantle grains as can be seen in Figures 
5(a) and 5(b). 
(b) V Hya 

V Hya was classified as C6,5 by Yamashita (1972). It 
is a variable star with overlapping periods: a period of 
about 530 days with amplitude of about 1.5 mg and a 
longer period of 6500 days with amplitude of 3.5 mg. This 
object is surrounded by an extended expanding molecu- 
lar envelope, resulting from extensive mass-loss. Mass-loss 
rate is in the range of 3.0 - 4.0 10 ~ 6 M /yr (Knapp 
& Morris 1985, Olofsson et al. 1990). Polarimetric mea- 
sures have been obtained by Johnson & Jones (1991) and 
more recently by Trammcll et al. (1994). Johnson & Jones 
(1991) have classified V Hya as a proto-planetary neb- 
ula and measured P(V) = 0.75 % ± 0.02% at 0=21° ± 
1°. Trammell et al. (1994) have observed this object in 
April 1992 and January 1993, and found that the polar- 
ization varied over this interval. The envelope properties 
found from molecular line observations by Knapp et al. 
(2000), like the fast molecular wind and the high mass 
loss rate, suggest that V Hya has entered its 'superwind' 
phase. However, its spectral type, period, colors, and lack 
of ionizing radiation indicate that this star is still on the 
AGB. Then, V Hya is believed to be in the latest phases 
of mass loss on the AGB. 

Our results show that the temperature of the cen- 
tral star is 2050K. Contrary to what is expected for this 
evolved phase, we have found a thin (r = 0.7) and exten- 
sive dust envelope (Ri=5.5R* and R 2 = 10000R,)- The 
size of the amorphous carbon mantle (niA.c.) is 950 A and 
the a-SiC core is 270 A. SiC/A.C. ratio is very small, 0.02. 
The best fit can be seen in Figures 5(c) and 5(d). We can 
say that a core/mantle grain model fits this star well, but 
a two homogeneous grain model cannot be discarded. 

SiC+: Stars 

We analyzed 7 out 11 SLMP objects classified as 
SiC+:, where the colon means a more uncertain classifica- 
tion. Almost all stars are well reproduced by a core/mantle 
grain model. The only exception is AX Cyg, which seems 
to need single particles of a-SiC, amorphous carbon and 
core/mantle grains simultaneously. The results obtained 
with the core/mantle grain code are similar to those found 
for SiC+ class stars. The temperatures of the central stars 
vary from 2400K to 2800K. Inner radii vary between 3 to 
4.9 R*. Carbon mantle sizes (mA.c.) vai T from 400A to 
1500 A with core (csic) values between 236A and 76lA. 
Optical depths are lower: 0.02 < r < 0.08. On the other 
hand SiC/A.C. abundance ratios vary from 0.10 to 0.27. 
(For these stars we could say that the core/mantle grain 
models are more adequate, even with such small optical 
depths, since the SiC/A.C. ratios are higher.) Figure 6 
shows an enlarged view of the two stars in this class, Y 
CVn and W CMa. We can see that the emission feature is 
shifted to longer wavelengths. This is a result of Mie's the- 
ory applied to core/mantle spheres. Suh (2000) has found 
a similar behavior in his models. 
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Fig. 5. This figure shows the best models for IRC+10216 
and VHya which belong to Red class. In Figure 5(a)-(d) 
we plotted models using core/mantle grain (solid line) 
and two homogeneous grain model (dashed line). Trian- 
gles represent the photometric data. In Figure 5(b) and 
5(d) we show an enlarged view of the 11.3 fim emission. 



(a) Y CVn 




7.2 7.8 8.4 9 9.610 12 

Wavelength (jim) 



Fig. 6. This figure shows best models for Y Cvn and W 
CMa which belong to SiC+: class. In Figure 6(a-b) we 
plotted models using core/mantle grain (solid line) and 
two homogeneous grain model (dashed line). Triangles 
represent the photometric data. 



4-2. Extreme Carbon stars 

Except for one case (R For), all SiC classes of the SLMP 
sample contains carbon stars which have thin dust en- 
velopes (r < 2.5). Our sample of extreme carbon stars 
contains objects that have optical depths varing between 
3.0 and 13. In order to make our analysis easier, we have 
decided to separate them according to this physical quan- 
tity: 3.0 < t < 5.2 (El group) and 7.0 < r < 13.0 (E2 
group). 

The stars belonging to our El group were better de- 
scribed by a two homogeneous grains model. Tempera- 
tures of the central stars are between 2200K and 2450K. 
Inner radii are about 6.5R, and outer radii 1000R* for all 
stars. Core (csic) and mantle (niA.c.) sizes are respectively 
about 200 A and 700 A in almost all cases. The abundance 
ratios SiC/A.C. vary from 0.03 and 0.07. These results can 
indicate that such stars are related to the SzCSLMP class. 

On the contrary, stars belonging to our E2 group were 
better fitted with the core/mantle grain code and are sim- 
ilar to the stars in the Red class. They are cooler (1900K < 
T e ff < 2200K) and present a dust envelope more extensive 
than those of the El group. Core and mantle sizes vary 
between 190 to 270 A and 700A to 1000A, respectively. 
The SiC/A.C. abundance ratios are also low, with values 
between 0.01 and 0.05. In the E2 group, three stars were 
better represented by taking into account a-SiC, amor- 
phous carbon and core/mantle grains simultaneously. We 
can speculate that these three stars represent a transition 
phase between SiC and Red classes. Figure 7 shows an en- 
larged view of AFGL 954 (El group) and AFGL 809 (E2 
group). 



5. Discussion and conclusions 

As commented in the previous section, the IR fluxes and 
the feature around 11.3 /zm present in most sources that 
we have analyzed, are likely to be better reproduced using 
two homogeneous grains (A.C. and a-SiC) simultaneously. 
This is true for SiC, SiC+ and our El stars. On the other 
hand, a few objects (those in Red and in our E2 class) are 
better described by a core/mantle grain model. However, 
the SiC ++ and Broad 1 classes cannot be reproduced 
with the existing optical constants. 

The SiC++ class presents both the 8-9^m emission 
and also a very prominent 14^tm absorption feature. The 
stars belonging to the Broad 1 class also show a 14^m ab- 
sorption. This absorption feature is weaker in Mira vari- 
ables than in SR variables. This may be explained by a 
stronger dust emission in Mira variables which fills the 
molecular absorption (see Yamamura et al. 1998). Regard- 
ing the origin of the 8-9/im emission feature, Aoki, Tsuji 
& Ohnaka (1999) have suggested that it may be a re- 
sult of molecular absorption at 7.5/^m and SiC emission 
at 11.3/im. This absorption could be due to HCN and/or 
C2H2 photospheric absorption bands. On the other hand, 
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Fig. 7. This figure shows the best models for AFGL 
954 and AFGL 809 which belong to El and E2 classes 
respectively. In Figure 7(a-b) we plotted models using 
core/mantle grain (solid line) and two homogeneous grain 
model (dashed line). Triangles represent the photometric 
data. 

the absorption feature at about 14/xm was attributed by 
the same authors to HCN and C2H2 absorption in the pho- 
tosphere or in the warm envelope close to the star. This 
absorption feature would be formed in the inner envelope 
where the mid-infrared radiation originates. On the other 
hand, Yamamura et al. (1998) have studied the 14/zm ab- 
sorption in the ISO SWS spectra of 11 carbon stars with 
mass-loss ranging from 10~ 8 to 10~ 4 M Q /yr. According to 
these authors, all stars clearly show an absorption band 
at about 13.7/zm due to C2H2 while the contribution from 
HCN molecules is small in this region. 

The stars of the SiC++ class have thin envelopes that 
can favor the misidentification of spectral features like 
the 8-9/nm one. On the other hand, this feature is also 
observed in some dust-enshrouded carbon stars, such as 
IRAS15194-5115, IRAS18239-0655 and IRAS18240+2326. 
Consequently, this feature could perhaps be produced by a 
solid component. Indeed, Goebel, Cheeseman & Gerbault 
(1995) have suggested a:C-H, as described by Dischler, 
Bubcnzcr & Koidl (1983). Unfortunately, they have not 
published the set of optical constants to this component. 
In order to try to reproduce the 8-9/zm emission, we have 
modeled SS Vir, considering optical constants of the amor- 
phous HAC as tabulated by Zubko (1996), but we have not 
been able to fit the emission. Bagnulo, Doyle & Andretta 
(1998) have shown that this feature can be reproduced by 
a dust shell composed of a mixture of SiC grains and sil- 



icate grains. This interpretation, however, does not seem 
to be consistent with the theory of dust formation. Be- 
fore discussing such results, let us comment now on some 
alternatives that might be investigated. 

On the basis of a simple model of circumstellar en- 
velopes, Kozasa et al. (1996) have proposed that the 
11.3/im feature could be attributed to small spherical 
core/mantle type grains (composed by a a-SiC core and 
a carbon mantle) in most cases. In their calculations they 
have used optical constants for SiC tabulated by Choyke 
& Palick (1985), which peak at about 10.7^im. In fact, 
using this set of constants, the emission produced by a 
core/mantle grain as proposed by them is about 11.3/im, 
and almost all sources of the SLMP sample could be fit- 
ted using a core/mantle grain code. This is no longer true 
when we use the constants proposed by Pegourie (1988), 
which peak at about 11.3/xm. In this case, the emission is 
shifted to longer wavelengths, too. We have computed our 
models considering the constants by Pegourie (1988) be- 
cause the full treatment consisting of the Kramcrs-Kronig 
analysis has been taken into account. 

Another possibility was raised by Speck et al. (1999): 
they have fitted some carbon stars using /3-SiC grains. Sil- 
icon carbide grains can be divided into two basic groups: 
a-SiC if the structure is one of the many hexagonal or 
rhombohcdral polytypes, and /3-SiC if the structure is cu- 
bic. /3-SiC feature occurs at about 0.4 /im shortwards of 
that of a-SiC. Their results were obtained without the 
KBr correction and they determined that a-SiC has an 
intense, broad band near 11.8/^m and /3-SiC peaks at 11.3 
to 11.4/zm. Silicon carbide grains found in meteorites have 
isotopic compositions that imply that most of these grains 
were formed around carbon stars. All studies to date of 
mctcoritic SiC grains have found them to be of the /3-type 
(Bernatowicz 1997). /3-SiC will transform into a-SiC above 
2100°C but the reverse process is thermodynamically un- 
likely The results obtained by Speck et al. (1999) show 
that there is an obvious predominance of the /3-SiC phase 
and that there is now no evidence for the a-SiC phase at 
all. Their sample contains SiC, SiC+:, Red and Extreme 
Carbon stars. However they do not solve the radiative 
transfer in these media. Moreover, we should expect some 
difference between "early" and "late" carbon stars with re- 
gard to dust grains. Unfortunately, the optical constants 
for /3-SiC were calculated in a short range of wavelengths, 
about 7 to 12 /zm. In this case, we need to adopt another 
set of optical constants at shorter wavelengths, where most 
of the stellar radiation is concentrated. With this assump- 
tion we cannot prove that /3-SiC grains are responsible for 
the 11.3 /im emission alone. 

SLMP have suggested the following carbon-rich dust 
sequence: SiC-f — ► SiC — > Red. The Red sources are signifi- 
cantly cooler on average than the SiC+ sources. Following 
our results, SiC+ stars have thinner envelopes than SiC 
stars. Temperatures of the central stars are very similar 
but there is a tendency to cooler temperatures in this se- 
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qucnce. Best models were obtained with two homogeneous 
grains for both SiC+ and SiC class. Red stars were best 
described with core/mantle grains. As mentioned before, 
based on our results we suggest that our sample of extreme 
carbon stars contains SiC and Red stars. In this sample, 
the thinner envelopes were best represented by two homo- 
geneous grains models (our El group) while thicker ones 
by core/mantle grain models (E2 group). The temperature 
of these stars are also cooler than in the SiC+ class. These 
results suggest that the sequence proposed by SLMP can 
be interpreted as an evolutionary scenario. Moreover, it 
seems reasonable to include our sample of extreme car- 
bon stars in such a scenario. In the beginning of the se- 
quence, grains are mainly composed of a-SiC and amor- 
phous carbon; with dust envelope evolution, carbon grains 
are coated in a-SiC ones. (Hence the emission is shifted 
to longer wavelengths). This phenomenon could perhaps 
explain the small quantities of SiC grains observed in the 
interstellar medium. 

Concerning SiC++ class stars, SLMP have proposed 
that they lie in a different evolutionary sequence, related 
to J-type carbon stars. In their sample, which contains 96 
objects, only nine are J-type stars. Two of them are classed 
in the SiC++ class. One of us (Lorenz-Martins 1996) have 
proposed that J-type carbon stars have an alternative evo- 
lutionary scenario which differs from that proposed for 
ordinary carbon stars. In fact, according to the results of 
the present paper, SiC++ stars have thicker (0.02 < t < 
0.20) envelopes than those of J-type stars (0.01 < t < 
0.05, see Lorenz-Martins, 1996). It seems that some cor- 
relation between SiC/A.C. ratios of both groups of stars 
also exists. Such results reinforce the suggestion by SLMP 
linking SiC++ and J-type carbon stars. 
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Table 3. Results of first sample 



Stars 


T cS (K) 


Ri (R*) 


R2 (R*) 


csiC (A) 


mA.c. (A) 


T 


SiC/A.C. 


Env. Class. 


00172+4425 


2200 


5.0 


1000 


453 


1000 


0.02 


0.10 


SiC++ 


01246-3248 


2400 


5.0 


1000 


500 


1200 


0.10 


0.08 


SiC++ 


02270-2619 


2500 


5.0 


5000 


407 


1050 


7.00 


0.06 


SiC 


03075+5742 


2400 


5.0 


800 


453 


1000 


0.01 


0.10 


N: 


03374+6229 


2650 


4.6 


1000 


500 


1200 


0.50 


0.08 


SiC+ 


04459+6804 


2700 


4.8 


1000 


500 


850 


0.03 


0.26 


SiC+: 


04483+2826 


2650 


3.0 


1000 


236 


400 


0.02 


0.26 


SiC+: 


04573-1452 


2250 


5.0 


1000 


294 


700 


0.60 


0.08 


SiC 


05028+0106 


2650 


5.2 


1000 


332 


700 


0.10 


0.12 


SiC+ 


05418-4628 


2400 


5.0 


1000 


500 


1200 


0.10 


0.08 


SiC++: 


05426+2040 


2600 


4.5 


1000 


554 


1000 


0.20 


0.20 


SiC 


05576+3940 


2200 


4.6 


1000 


403 


1000 


1.20 


0.07 


Br2 


06225+1445 


2700 


4.9 


1000 


400 


850 


0.03 


0.12 


Br2 


06331+3829 


2550 


3.7 


1000 


340 


700 


0.10 


0.13 


SiC 


06529+0626 


2200 


4.6 


1000 


282 


700 


1.00 


0.07 


SiC 


07045-0728 


2400 


6.5 


1000 


270 


500 


0.10 


0.18 


SiC+ 


07057-1150 


2650 


3.0 


1000 


440 


740 


0.04 


0.27 


SiC+: 


07065-7256 


2400 


4.6 


1000 


282 


700 


2.20 


0.07 


SiC 


08538+2002 


2400 


5.0 


1000 


400 


1200 


0.20 


0.04 


SiC++ 


09452+1330 


2100 


5.5 


8000 


110 


500 


10.0 


0.01 


Red 


10329-3918 


2700 


4.5 


1000 


390 


850 


0.04 


0.11 


SiC+: 


10350-1307 


2700 


4.5 


1000 


468 


850 


0.03 


0.20 


SiC 


10491-2059 


2050 


5.6 


10000 


270 


950 


0.70 


0.02 


Red 


12226+0102 


2700 


5.0 


1000 


487 


1200 


0.35 


0.07 


Brl 


12427+4542 


2700 


4.9 


1000 


317 


700 


0.05 


0.10 


SiC+: 


12447+0425 


2200 


4.3 


1000 


307 


700 


2.50 


0.09 


SiC 


12544+6615 


2650 


3.7 


1000 


155 


400 


0.04 


0.06 


SiC+ 


15094-6953 


2650 


5.3 


1000 


416 


700 


0.03 


0.26 


SiC+ 


18306+3657 


2200 


5.0 


1000 


317 


700 


0.03 


0.10 


SiC++ 


19017-0545 


2550 


4.9 


1000 


315 


800 


0.10 


0.07 


SiC+ 


19555+4407 


2400 


4.0 


1000 


761 


1500 


0.08 


0.15 


SiC+: 


21032-0024 


2200 


4.5 


1000 


294 


700 


2.50 


0.08 


SiC 


21399+3516 


2800 


4.0 


10000 


500 


1000 


0.04 


0.14 


SiC+: 


23587+6004 


2500 


3.0 


800 


554 


1000 


0.01 


0.20 


Brl 



Table 4. results of the Sample with evolved stars 



IRAS 


T eB (K) 


Ri (R*) Ri (R*) 


CSiC (A) 


mA.c. (A) r 


SiC/AC 


Env. Class 


05377+1346 


2350 


3.0 


2000 


282 


700 


3.0 


0.07 


El 


05405+3240 


2200 


7.1 


10000 


215 


700 


10.6 


0.03 


E2 


06012+0726 


2000 


4.5 


4000 


189 


900 


12.5 


0.01 


E2 


06291+4319 


2300 


6.5 


1000 


226 


700 


3.5 


0.04 


El 


06342+0328 


2500 


4.0 


7000 


256 


700 


7.0 


0.05 


E2 


07098-2112 


2450 


6.5 


1000 


236 


700 


4.0 


0.04 


El 


08088-3243 


2300 


6.4 


1000 


215 


700 


4.0 


0.03 


El 


15082-4808 


2050 


4.6 


1000 


273 


1000 


9.4 


0.02 


E2 


19594+4047 


2000 


5.0 


7000 


232 


800 


13.0 


0.025 


E2 


20570+2714 


2200 


8.0 


1000 


156 


500 


5.2 


0.03 


El 


23257+1038 


1900 


7.0 


1000 


254 


700 


13.0 


0.05 


E2 



